Spin dynamics in tunneling decay of a metastable state 
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We analyze spin dynamics in the tunneling decay of a metastable localized state in the presence 
of spin-orbit coupling. We find that the spin polarization at short time scales is affected by the 
initial state while at long time scales both the probability- and the spin density exhibit diffraction- 
in-time phenomenon. We find that in addition to the tunneling time the tunneling in general 
can be characterized by a new parameter, the tunneling length. Although the tunneling length is 
independent on the spin-orbit coupling, it can be accessed by the spin rotation measurement. 
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I. INTRODUCTION 

Spin-orbit (SO) coupling, interaction of particle spin or 
pseudospin with the orbital motion, provides an efficient 
way to control and manipulate spin, charge, and mass 
transport. Two different kinds of SO coupled systems at- 
tract a great deal of interest due to the known and yet un- 
explored variety of phenomena they can demonstrate and 
due to their possible applications in quantum technolo- 
gies. One class of systems, investigated for more than five 
decades by now, is semiconductors and semiconductor- 
based nanostructures [H-Ql- The other class, lavished at- 
tention only recently, is cold atoms and Bose-Einstein 
condensates, where by engineering external optical fields 
one can couple orbital motion to the pseudospin degree of 
freedom |8l-[l3| and cause Dresselhaus and Rashba types 
of SO coupling similar to that in solids. 

In quantum systems of interest the tunneling plays 
an important role and either completely determines or 
strongly influences the particle dynamics. At certain con- 
ditions the tunneling rate depends on the spin of the 
particle 0, S [H, Ej]- The understanding of the tun- 
neling is the key for the understanding of the transient 
processes in a broad variety of systems, including, e.g. 
charge transport in molecular nanostructures [T5| . One 
of the key issues in the tunneling theory is the evalua- 
tion of the time spent by the particle in the classically 
forbidden regions. Similar to the scattering problem of 
propagating wave packets, the question of the tunnel- 
ing time for the decay of a metastable system [T6l - [20j 
may also be posed. The spin-dependent effects based on 
the Larmor clock concept for potential barriers [111, , 
Hartman effect in graphene [23|, and effective exchange 
fields in semiconductors can provide a measure of 
this time. In classically forbidden regions there is not 
only precession but also a rotation of the magnetic mo- 
ment into the magnetic field direction and the time of the 
interaction with the barrier is closely related to this rota- 
tion [2^, [2^ . The effects of SO coupling on the tunneling 
through semiconductor quantum-well structures with a 
lateral potential barrier provide another tool to 

utilize electron spin modifying the charge transport in 



nonmagnetic systems. 

Here we investigate the spin-dependent tunneling of 
a state initially localized in the potential at short and 
long compared to the state lifetime time scales [2^ . This 
paper is organized as follows. In Sec. II we introduce 
the model Hamiltonian and formulate the physical ob- 
servables of interest. We use the SU(2) spin rotation to 
gauge out spin-orbit coupling and restore its effects in the 
calculation of the observables by inverse transformation. 
In Sec. Ill we study the dynamics at short time scale, 
concentrating on spin oscillations, decay, and escape from 
the localizing potential. In Sec. IV we study long-term 
dynamics in the far-field zone. The diffraction in time in 
spin density is observed at a long distance, where SO cou- 
pling forms a precursor in the propagating density. We 
show that the tunneling cair be characterized not only by 
time, but also by a certain length parameter, which can 
be accessed by detecting the spin precession due to the 
SO coupling. Conclusions summarize the results. 



II. MODEL FOR SPIN-DEPENDENT 
TUNNELING 

As a model we consider shown in FiglT]time-dependent 
potential U{x,t), infinite at x < 0, with the time- 
dependence: 



(1) 



U{x,t) = Uiix) (i<0), U{x,t)^U2ix) {t>0). 

At i < the potential holds bound states with wave func- 
tions fj (x) and energy Ej ; it changes at t > from a step 
to a barrier to allow the tunneling. Such time-dependent 
potential can be produced by a recently developed tech- 
nique [13 where a moving laser beam "paints" a broad 
variety of coordinate- and time-dependent potentials. 

With SO coupling taken into account, the total one- 
dimensional Hamiltonian is 
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A, 



C/(.)-f^, (2) 



with k 

Pso 



= —id/dx, the vector potential ~ Pso^x, 
ma/h, where a is the Dresselhaus SO coupling 
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FIG. 1. Time-dependent potential U{x): Ui{x) at t < and 
U2{x) when t > 0. Ui{x) is a step in the positive halfplane, 
that is zero from a- = to a; = ai and Uo from x = ai to 
infinity, while U2{x) is a barrier of the height Uo extended 
from X = ai to a; = 02. The potential is always infinite in the 
negative halfplane. The barrier width is defined as d = 02— ai. 



constant, and is the corresponding Pauli matrix. The 
spatial scale of spin precession is characterized by the 
length 2^ = h/pso- The typical values of f for differ- 
ent systems vi^itli SO coupling can be of interest. For (i) 
electrons in semiconductor GaAs nanostructures [l| with 
m « 6 X 10~^^ g, a/h ^ 0.5 x 10^ cm/s, (ii) for ^Li atoms 
[il with m ss 10-23 g, a/h ^ 10 cm/s, and (iii) for ^^Rb 
atoms [iJl with to w 1.5 x 10-^^ g, a/h ^ 0.3 cm/s, we 
obtain ^ ^ 10"^ cm. It is interesting to mention that 
although masses and coupling constants for electrons in 
nanostructures and cold atoms are very different, the re- 
sulting precession length is on the same order of magni- 
tude for all these systems. Since typical localization scale 
of electrons in nanostructures oi is on the order of 10"^ 
cm, they are in the weak-coupling regime with ^/ai 3> 1. 
Cold atoms, however, are localized on the scale of the or- 
der visible light wavelength, that is of IQ-* cm, and the 
strong SO coupling regime with ^ < ai can be achieved 
there. The characteristic timescale corresponding to the 
particle motion inside the potential, ma\/h, is on the or- 
der of 0.1-1 ps for electrons in nanostructures and 0.1-1 
ms for cold atoms. 

From now on we use the system of units with ft = 1 
and particle effective mass to = 1. The wave function 
corresponding to H is ij}{x^ t) with the initial state set as 



(3) 



and the spinor x = [IjO] corresponds to the z-axis ori- 
entation. Such a function can be produced by applying 
a magnetic field along the z-axis for an electron or by 
a special design of optical field in the case of a neutral 
bosonic atom [T^]. At t > 0, neither the orbital wave 
function tpj is the eigenstate, nor x is the eigenspinor 
of Hamiltonian H . The initial state begins to evolve at 
t > with spin precession and the probability to find the 
particle inside the potential decreases simultaneously. 

We use a gauge transformation H = SHS^^ with a 
unitary matrix S = expiia xX / 2^/) to gauge away the SO 
coupling. The gauge transformation shifts by Ox/^^ 



and turns SO coupling into a constant, with the wave 
function evolving in time transformed as: 



= Sip{x,t) 
with the upper and lower components 



■ipi{x,t) 
4>-i{x,t) 



'ip^{x,t) = / G^(fc)0fe(a;)exp 
Jo 



(4) 



dk, (5) 



and the coefficients are 
G.(fc) = 



^^}a■ix, 0)(f)kix)dx. 



(6) 



Here, (pkix) is the eigenstate of the Hamiltonian corre- 
sponding to U2{x): 

!C{k) siw{kx), {Q < X < ai) 

^(fc)e^K,. ^ F(fc)e«'=^ {ai<x< 02) (7) 
y277rsin[fca; + ^(fc)], (x > 02) 

normalized as (0^/ — 5{k~k'). The coefficients C(fc), 
D{k), F{k) and the phase 9{k) satisfy the boundary con- 
ditions of the potential U2{x), and Kk = V^Uq — k"^. 
In the tunneling regime k < \/2Uo, while in the prop- 
agating regime k > ^/2Uo, and ink is substituted by 
q = Vk'^ - 2f/o. 

The initial wave function corresponding to Hamilto- 
nian H can be expressed as 



il){x,0) = ip-i{x) 
The coefficients Gcr(fc) become 



cos(x/2^) 
i sin(x/2^) 



(8) 

(9) 
(10) 



In what follows, we investigate the particle motion by 
calculating the physical observables such as probability 
density, spin density, and spin polarization defined as: 



Gi(fc) = / coa{x/2^)ipj{x)(f>k{x)dx, 
Jo 

POO 

G-i{k) = i sm{x/2S^)ipj{x)(f>k{x)dx. 



p{x,t) = 'il)''{x,t)ip{x,t), 
ai{x,t) = ■4){x,t)''d-^ip{x,t), 
d-i{x,t) 



p{x,t) 



(11) 
(12) 

(13) 



respectively. As the spin rotates around the effective 
SO coupling field along the a;-direction, the integral of 
CTa;(a;, t)-component over the .t > half axis is conserved. 
For this reason, we investigate spin density and spin po- 
larization in the more informative y— component. 
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III. SPIN DYNAMICS AT SHORT TIMES 

To address the short-term dynamics on the time less 
than the hfetime of the initial bound state [2^ , we study 
two relevant quantities. First quantity is 17^(02, i), the 
spin density at the exit of the barrier 02- The other one, 
Py™' (t) , defined as 



il}\x, t)ayip{x, t)dx 
t/j^(.T, t)dx 



(14) 



being the spin polarization y component in the poten- 
tial, represents the integrated spin dynamics inside the 
potential. The denominator of Eg. p^ is the probability 
to find the particle inside the potential, which decays to 
1/e of its initial value at the lifetime of the metastable 
state. 

We consider the evolution of three initial orbital states, 
as shown in FigH] for the spin component at the edge. 
We choose Uo = 16 so that there are two bound states, 
the ground state (y9o and the first excited state (/Si, and 
use dimensionlcss ai = 1. Here, the initial state has 
a strong impact on the time-dependent (jy{a2,t). It is 
demonstrated in Figl2] that it takes very short time for 
ay{a2,t) to develop into the minimum from zero, irre- 
spective of the initial state. This behavior is similar to 
the fast development of a. plateau in the outgoing fiux af- 
ter the potential change [29[. However, spin density with 
ifi decays faster than that with ipQ. The linear combi- 
nation of these two bound states (pcom — i'^o + '^i)/V2 
shows strong oscillations due to the interference between 
them. As both spin density and probability density in 
the potential decay with time, spin polarization tends to 
a constant at large time. 




FIG. 2. (Color online) Comparison of time evolution of 
ay{a2,t) at the right edge of the barrier 02 with different ini- 
tial states, the ground state ifio (solid), the first excited state 
ifi (dashed) and the linear combination 

(dot-dashed), provided by the barrier d = 0.4, Uo = 16, and 
^ = 0.5. Negative ay{a2,t) is related to the direction of spin 
precession determined by the sign of coupling constant a. 

The tunneling resulted from various initial states has 
a strong effect on spin polarization, especially in short- 



time scales, as shown in FiglS) Spin polarization inside 
the potential oscillates between boundaries. The oscilla- 
tion rate is fast and determined by the energy difference 
of the initially bound states. With the time, the con- 
tribution of the excited state rapidly decays due to the 
fast tunneling, and the spin remains in the state achieved 
by the time of the decay of the upper bound state. The 
amplitude of resonances for ipi is larger than ipo, as the 
former possesses larger momentum. As for the linear 
combination of these two bound states, the spin polar- 
ization is greatly enhanced because of interferences. On 
the other hand, strength of SO coupling strongly infiu- 
ences the spin polarization in the potential as illustrated 
by making comparisons for ^ = 0.5 and ^ = 1 parameters. 
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FIG. 3. (Color online) Time evolution oi Py (t) in the poten- 
tial from to ai, caused by SO coupling with ^ — 0.5 (solid) 
and f = 1 (dashed), with different initial coordinate states: 
(a) ifio, (b) ifii, and (c) i^com = {<po + fi) / V^, provided by the 
same barrier as in Fig[2] 

Having illustrated the evolution of spin density at the 
boundary, ay{a2,t), and averaged spin polarization, we 
can address the details of the spin distribution inside the 
potential. As shown in FigHl this distribution oscillates 
due to the interferences between the ground state and 
other eigenstates. The amplitude of oscillations becomes 
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weaker with time as the state decays. 
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FIG. 4. (Color online) Spin density {y component) in the 
potential vs time and position for ^ = 0.5 and i]}{x,t = 0) = 
Xfo- The parameters of the barrier are the same as in Fig[2] 



IV. SPIN DYNAMICS IN THE FAR-FIELD 
ZONE 



Here we investigate the long-term spin dynamics, by 
considering the probability density p and y component 
of spin density Uy (X, t) which can be detected at a given 
long distance X ^ 1. The density evolution for p{X,t) 
a.t X ~ IOtt is shown in FigE] for different SO coupling 
strength. Under different ^, the time evolutions of prob- 
ability density show a strong sharp peak followed by os- 
cillations due to the diffraction-in-time process 




However, unlike the case without SO coupling [2£ 
for strong couplings (smaU ^), there exist some oscilla- 
tions before the sharp peak, where the interference of two 
velocities of spin up and spin down components are re- 
markable. A very interesting feature is the precursor of 
the main peak corresponding to the opposite spin. The 
precursor becomes stronger with the increase in the SO 
coupling. The time dependence of ay{X,t) ed, X = IOtt 
is shown in FiglHl Different values of ^ result in differ- 
ent spin rotation angles and corresponding spin density, 
provided that ay{x,i) is detected at the same position. 
Spin precession in the (y, z)-plane is described by the 



classical precession angle, /S, 



[so] 



fJi, where the preces- 



sion frequency is 2afco = fco/'? f^d the wave packet 
is detected at the position X at the time instant X/v, 
the velocity w = fep. Therefore, the rotation angle is 
/3g^°' = X/^, independent on time t, provided that the 
wave packet is considered classically. Ideally, for a free 
particle propagating from the origin x{t = 0) = 0, spin 

polarizations components are Px = 0, Py = — sin/Jp^"', 

Pz = cos/3q''°'. As a result, at X = mr^, the classical spin 
polarization should be Py = 0. However, in the tunneling 
problem we consider, the value of Py is not zero even at 
X = niT^. 




FIG. 5. (Color online) (a) Time-dependent density at the 
observation point X — IOtt for system parameter Uo = 16, d = 
0.4, and ^ = 0.5. Inset shows the precursor dynamics on a 
shorter time scale. Dashed line corresponds to Uo = 16, d = 
0.4, and ^ = 1. (b) Solid line corresponds to Uo — 16, d = 0.4, 
and ^ = 0.5, dot-dashed line corresponds to {/o = 8, d = 0.4, 
and ^ — 0.5. Inset shows the precursor dynamics on a shorter 
time scale. 
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FIG. 6. Spin density in the observation point X = IOtt 
for Uo = 16, d = 0.4, %p{x,t = 0) = xvo and C = 0.5. The 
precursor is made by the contribution of the opposite spin. 
Other parameters are the same as in Fig[2] 



For the weak SO coupling ^ » 1, where <C 1, 

the spin polarization can be calculated as: 



Py Ri — sm — + 2 cos — ^ — 

C ^ ■01 



(15) 



Second term in Eg. p^ can be viewed as a result of a 
correction to the particle displacement in the form py = 



5 
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FIG. 7. (Color online) Parameter of additional displacement 
Sx as the function of the barrier transparency. Filled circles 
correspond to Uo — 16 and variable width d, the squares 
correspond to d = 0.4 with variable Uo, and ^ = 10. 



— sin(.T/^ + 5x/£,), where, according to Eas.(P)). pIH) : 

poo 

/ G_i(fc)exp(-ifc2f/2)(/)fe(a;)dfc 
Sx = -2 lim ^^^^^ ^ (16) 

Jo 

This deviation of Py from zero at X ~ nvr^ demonstrates 
that the real precession angle is ~ /3q°' + 

with the correction 6/3^^°^ can be viewed as a result of ad- 
ditional displacement of the particle dx with Sx = ^5/3 [''"I . 
From Eg. p^ . we can see that Sx is independent of the 
SO coupling, while the parameters of the barrier play an 
important role, as the cigenfunctions 4>k{x) strongly de- 
pend on them. Numerical results also show that Sx is 
independent of coordinate and time if measured at dis- 
tance X ^ a2 at times t > X/v. The tunneling length Sx 
is presented in FiglT] for two different barriers. This fig- 
ure shows a clear crossover from the classical (transparent 
barrier) to the tunneling (opaque barrier) regimes, where 
Sx shows a saturation, and the additional displacement 



is universal. The saturation of the tunneling length for 
opaque barriers seems similar to the Hartman effect (35j 
on the group delay in a scattering case [s^. However, 
after making comparisons of St = Sx/v and calculated 
group delay for a single barrier with parameters presented 
in FiglZl we cannot draw the unambiguous conclusion on 
the relevance of these two quantities. 



V. CONCLUSIONS 

In conclusion, we investigated the spin dynamics of 
a tunneling particle initially localized in a potential in 
the presence of SO coupling. It is shown that at short- 
time scales initial states play an important role in spin 
polarization while the spin density and probability den- 
sity possess diffraction-in-time phenomenon at long-time 
scales. We showed that in addition to the time, tunnel- 
ing can be characterized by a characteristic length. We 
use the rotation angle to identify the tunneling length in 
the presence of weak coupling. The tunneling length de- 
pends on the barrier parameters, being independent on 
the strength of the spin-orbit coupling. These effects can 
be observed in experiments with cold atoms, where the 
SO coupling is strong enough to cause spin precession on 
a relatively short spatial scale. 



ACKNOWLEDGEMENT 

This work of FYS was supported by the MCINN of 
Spain Grant FIS2009-12773-C02-01, by "Grupos Consol- 
idados UPV/EHU del Gobierno Vasco" Grant IT-472-10, 
and by the UPV/EHU under program UFI 11/55. Y. B. 
acknowledges financial support from the Basque Govern- 
ment (Grant No. BFT2010-255). We are grateful to J.G. 
Muga and D. Sokolovski for valuable discussions. 



[1] 1. Zutic, J. Fabian and S. Das Sarma, Rev. Mod. Phys. 

76, 323 (2004). 
[2] S. Amasha, K. MacLean, I. P. Radu, D. M. Zumbiihl, 

M. A. Kastner, M. P. Hanson, and A. C. Gossard, Phys. 

Rev. B 78, 041306(R) (2008). 
[3] P. Stano and P. Jacquod, Phys. Rev. B 82, 125309 

(2010). 

[4] P. Stano and J. Fabian, Phys. Rev. B 72, 155410 (2005). 
[5] R. Romo and S. E. UUoa, Phys. Rev. B 72, 121305(R) 
(2005). 

[6] 1. Malajovich, J. M. Kikkawa, and D. D. Awschalom, J.J. 

Berry and N. Samarth, Phys. Rev. Lett. 84, 1015 (2000). 
[7] D. D. Awschalom and M. E. Flatte, Nat. Phys. 3, 153 

(2007). 

[8] T. D. Stanescu, B. Anderson, and V. Galitski, Phys. Rev. 
A 78, 023616 (2008). 



[9] X.-J. Liu, M. F. Borunda, X. Liu, and J. Sinova, Phys. 

Rev. Lett. 102, 046402 (2009). 
[10] G. Wang, G. Gao, C.-M. Jian, and H. Zhai, Phys. Rev. 

Lett. 105, 160403 (2010). 
[11] Y.-J. Lin, K. Jimenez- Garcia, and I. B. Spielman, Nature 

471, 83 (2011). 
[12] D. L. Campbell, G. Juzeliiinas, and I. B. Spielman, Phys. 

Rev. A 84, 025602 (2011) 
[13] M. M. Glazov, P. S. Alekseev, M. A. Odnoblyudov, V. M. 

Chistyakov, S. A. Tarasenko, and I. N. Yassievich, Phys. 

Rev. B 71, 155313 (2005). 
[14] D. Khomitsky and E. Ya. Sherman, EPL 90, 27010 

(2010). 

[15] X. Zheng, F. Wang, C.-Y. Yam, Y. Mo, and G.H. Chen, 
Phys. Rev. B 75, 195127 (2007), Y. Wang, C.-Y. Yam, 
G.H. Chen, T. Frauenheim, and T.A. Niehaus, Chem. 
Phys. 391, 69 (2011), G.-Y. Yam, X. Zheng, G.H. Chen, 



Y. Wang, Th. Frauenheim, and T.A. Niehaus, Phys. Rev. 
B 83, 245448 (2011). 
[16] A. Marchewka and E. Granot, Phys. Rev. A 79, 012106 

(2009) . 

[17] G. Kalbermann, Phys. Rev. C 79, 024613 (2009); 77, 
041601 (2008). 

[18] N. G. Kelkar, H.M. Castaneda, and M. Nowakowski, Eu- 

rophys. Lett. 85, 20006 (2009). 
[19] S. Cordero, G. Garcia-Calderon, R. Romo, and J. Villav- 

icencio, Phys. Rev. A 84, 042118 (2011). 
[20] M. Pons, D. Sokolovski, and A. del Campo, Phys. Rev. 

A 85, 022107 (2012). 
[21] A. I. Baz', Yad. Fiz. 5, 229 (1967). 
[22] V. F. Rybachenko, Yad. Fiz. 5, 895 (1967). 
[23] R. A. Sepkhanov, M. V. Medvedyeva, and C. W. J. 

Beenakker, Phys. Rev. B 80, 245433 (2009). 
[24] B. Huang and I. Appelbaum, Phys. Rev. B 82, 241202(R) 

(2010) . 

[25] M. Biittiker and R. Landauer, Phys. Rev. Lett. 49, 1739 
(1982). 

[26] M. Biittiker, Phys. Rev. B 27, 6178 (1983). 
[27] V. A. Sablikov and Y. Ya. Tkach, Phys. Rev. B 76, 
245321 (2007). 



6 



[28] V. I. Perel', S. A. Tarasenko, I. N. Yassievich, S. D. 
Ganichev, V. V. Bel'kov, and W. Prettl, Phys. Rev. B 
67, 201304 (2003). 

[29] Y. Ban, E. Ya. Sherman, J. G. Muga, and M. Biittiker, 
Pliys. Rev. A 82, 062121 (2010). 

[30] K. Henderson, C. Ryu, C. MacCormick and M. G. 
Boshie, New J. Phys. 11, 043030 (2009). 

[31] M. Moshinsky, Phys. Rev. 88, 625 (1952). The result 
of this paper can be understood in terms of diffraction 
of perfectly coherent propagating waves. For decay of a 
metastable state, the diffraction can be seen as produced 
by partially coherent waves: M. Born and E. Wolf, Prin- 
ciples of Optics, Chapter X, 1999 (7th Ed.), see also in 
Ref. [3|. 

[32] A. Steane, P. Szriftgiser, P. DesbioUes, and J. Dalibard, 

Phys. Rev. Lett. 74, 4972 (1995). 
[33] A. del Campo, G. Garcia-Calderon, and J.G. Muga, 

Phys. Rep. 476, 1 (2009). 
[34] E. Torrontegui, J. Muiioz, Y. Ban, and J. G. Muga, Phys. 

Rev. A 83, 043608 (2011). 
[35] T. E. Hartman, J. Appl. Phys. 33, 3427 (1962). 
[36] E. P. Wigner, Phys. Rev. 98, 145 (1955). 



